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Abstract 

We study the ground state of the two-dimensional Anderson-Hubbard model 
using a quantum real space renormalization group method. We obtain the 
phase diagram near half filling. The system is always insulating with disorder. 
At half filling, the system undergoes a transition from a gapless (Anderson) 
insulator to an incompressible (Mott-Hubbard) insulator as the interaction 
U reaches a critical value U c . Away from half filling, the insulating phase 
is always gapless and is found to be controlled by the Anderson fixed point 
at half filling. This result is similar to that obtained in the corresponding 
one dimensional system and suggests strongly the importance of the electron- 
electron correlation in this gapless insulating phase. 
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I. INTRODUCTION 



When disorder is introduced into a physical system, it will result in the localization 
of some single particle states Q. For a noninteracting electron system, if the states at 
the fermi surface are localized, the system is a so-called 'Anderson insulator'. In reality, the 
interaction between electrons always exists, and such a single-particle picture may not apply. 
The understanding of the interplay between disorder and interaction has been an important 
issue in condensed matter physics H. This issue has become more interesting since the 
discovery of the high temperature superconductors O. On the one hand, it is commonly 
believed that these materials are strongly two dimensional (2D) in character, and that 
electron-electron correlations are important and are responsible for many of their unusual 
physical properties; on the other hand, it is also clear that disorder, which manifests itself 
as (e.g.,) oxygen vacancies, is inevitably present in these materials. It is thus of interest to 
investigate the effect of disorder in highly correlated systems. In this paper, we pursue such 
a study using a real space renormalization group (RG) approach. The model we consider is 
the two dimensional Anderson-Hubbard model defined by the Hamiltonian 

H = J2 (tijctcj* + H.C.) + J2( W i - f i ) n icr + U J2 n H n H- C 1 - 1 ) 

<i,j> icr i 

a 

Here c i(T (ci a ) is the creation (annihilation) operator for a spin-a electron on site i, ty is the 
nearest-neighbor intersite hopping energy and U (> 0) gives the on-site Coulomb repulsion 
energy. The chemical potential is given by fi, and Wi is a random site potential which 
has an independent gaussian distribution with zero mean and width W, i.e., W{ = and 
WiWj = W 2 5ij (where the overbar indicates random average). We shall consider only the 
square lattice case. Without interaction (U = 0) but with disorder, this is the Anderson 
model of localization which has been the prototype for studying the effect of disorder 
in electron systems. With interaction but without disorder (Wi = 0), this is the Hubbard 
model 0, which is believed to be one of the simplest theoretical models which possesses 
the essential physics of correlation ||], and which has been the focus of intense theoretical 
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investigation since the discovery of high T c superconductors. Thus, the Anderson- Hubbard 
model is a natural starting point for the investigation of the combined effects of disorder 
and interaction in electron systems. Throughout the paper our discussion will be restricted 
to ground state properties. 

For the Anderson model in 2D, the consensus is that all of its eigenstates are localized, 
and hence that it describes a gapless insulating state ||. For the Hubbard model at half 
filling, on the other hand, it is commonly accepted that the Coulomb repulsion U gives rise to 
insulating behavior, with long-ranged antiferromagnetic ordering in the ground state. Thus, 
in contrast to the usual 'band insulator' where the insulating phase is due to the filling of 
electron bands in the solid, and different from the Anderson insulator where the vanishing 
conductivity is caused by the localization of the single particle states at the fermi surface, 
the insulating state in the Hubbard model at half filling is a result of electron-electron 
correlations, hence a 'correlated insulator'. Away from but close to half filling, the Hubbard 
model describes a highly correlated system whose exact properties we still know little about 
despite intense studies during the past few years. Candidates for the possible ground states 
can be, e.g., phase separation |3|, a highly correlated metal 0, or a superconductor 0,0- 
When both interaction and disorder are considered, one expects that disorder breaks the 
translation and other lattice symmetries and possibly weakens effect of the correlations; 
on the other hand, strong correlations may render the standard single particle picture of 
Anderson localization meaningless. As a first step towards understanding this complicated 
issue, we wish to identify the phase diagram of the 2D Anderson-Hubbard model. 

Previously, the Landau Fermi-liquid idea has been employed to describe systems of 
weakly interacting electrons with (weak) disorder 0. The validity of approaches along 
this line is questionable in the present situation because the noninteracting system is non- 
metallic, and because the effect of the interaction is presumably nonperturbative. The real 
space RG approach ||- j[3|, on the other hand, is a nonperturbative method which allows 
one to treat disorder and interaction of any strength on the same footing. It has however 
the disadvantage of being an uncontrolled approximation, so that its implementation and 
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interpretation should be taken with extra caution. 

The real space RG scheme adopted in this paper is a generalization of the works of 
Hirsch []8| and Ma |9j . This method allows one to study the compressibility of the system 
by investigating the renormalization of chemical potential and the corresponding flow of 
density. This RG scheme has been previously employed to study the U = oo Anderson- 
Hubbard model for spinless bosons fll|Jl2|| . While the quantitative results, such as the 
critical exponents of the superfluid-Bose-glass phase transition, are still the subject of some 



controversy |L4]], this method does provide the correct qualitative physical picture. For 
instance, it shows that the superfluid phase is unstable against any amount of disorder in 
the ID U = oo Anderson-Hubbard model, in agreements with the exact result ||15|| : in 2D 



and 3D, it shows a transition from the superfluid phase to a disordered (Bose-glass) phase at 
some critical amount of disorder, as indicated by other theoretical approaches |14| . Thus we 



have reasons to believe that the real space RG approach can also give us useful information 
concerning the fermion Anderson-Hubbard model ( |1 . 1| ) . In addition, the validity of the real 
space RG scheme for the present case may be tested in the noninteracting system; for this 
case our method gives the result that the metallic phase is unstable against any amount of 
disorder (see below), consistent with the now accepted theoretical results @.|nj . 

One may wonder why we are in a position to investigate the Anderson-Hubbard model 
when there is still not a good understanding of even the pure system. Our response to this 
question is that, as discussed previously |11| , |12| , and as will be emphasized in the section to 



follow, disorder is in fact an advantage for our investigation. By sampling a large ensemble 
of random configurations of the potential {H^}, the problem of losing long range quantum 
correlations due to breaking the system into blocks in the real space RG may be partially 
compensated. Also, the disorder averaging allows us to treat the (average) particle n as 
a continuous variable — which is not possible in the absence of disorder. Finally, even for 
disordered systems there already exist some known cases where one can test the method, 
as discussed above. These considerations give us some confidence that our real space RG 
approach to the disordered problem is a suitable choice, at least for our present purpose of 



investigating the phase diagram. 

The rest of the paper is organized as follows: the real space RG method is described 
in some detail in the next section; in Section III, we present our results; and we offer a 
summary and discussion in Section IV. 



II. METHOD 



Our real space RG method is similar to that of Refs. || and JTT|. This real space RG 



is implemented numerically on a finite lattice. The random field {Wi} is obtained through 
a gaussian random number generator. Briefly, the RG procedure can be described in the 
following five steps: (i) Divide the lattice into blocks of size n s . (ii) Compute the block 
fermion operators which are defined in terms of four eigenstates of the block Hamiltonian. 
Each block is characterized by an effective on-site potential and an on-site repulsion between 
the 'block particles'. (Hi) Determine the hopping parameters for the block particles from 
the inter-block couplings between the site variables. However, since the block parameters 
arise from the random {W^}, and so have in general a different distribution from the original 
one, we need to (iv) repeat the above procedures [(ii) and (Hi)] for a large random ensemble 
to determine the distribution of the block parameters. We shall limit ourselves to tracking 
only the first two moments of each distribution (see below). At this stage, the Hamiltonian 
is mapped back onto its original form with renormalized parameters: 

H = (WcLc/3a + H.c.) + J2(W a - fl)n aa + ™^n al + constant (2.1) 

<a/3> aa a 

where a and (3 are block indices. Finally, (v) we the iterate above sequence to find the flow 
and fixed point (s) of the RG. 

Now we elaborate each of the above steps: 

(i). The blocks we used in the present work are shown in Fig. 1. Each of them is chosen 
to have an odd number of sites, to allow us to correctly treat the physics of half filling JT7 . 



Even for such small blocks the numerical diagonalization is non-trivial, due to the large 



Hilbert space and the necessity of sampling a large number of random configurations. For 
the 3x3 square lattice at half filling, for instance, the dimension of the Hilbert space is 
15876 x 15876. (Here the lattice symmetry is destroyed by disorder and thus can not be 
used to reduce the size of the Hilbert space.) Both types of blocks shown in Fig. 1 have 
been tested, and we find that they give the same qualitative physics. Thus for our present 
purpose (exploring the phase diagram) we may focus on the 'star' block (Fig. la) which is 
computationally more convenient. 

(ii). Each (microscopic) site can have one of four possible states: the no-electron state 
|0), the up-spin electron state the down-spin electron state ||), and the two-electron state 
III). The energy for the no-electron state is denoted and the two-electron state energy 
E( 2 \ The up- and down-spin electron state energies are degenerate and denoted by 
For each block, we find the exact ground state and ground-state energy for the Hamiltonian 
(Eqn. |0|) for every possible odd number of particles, restricted to the subspace with S = \ 
(the two S g — ±| states are degenerate). The lowest energy among these ground state 
energies then gives us E( l > . Letting N a be the number of particles which gives 

E^ in block 

a, and | | a ) the corresponding ground state, we then take the N a — 1 ground state (from the 
subspace S = 0) as IOq.) (with energy E^), and the N a + 1 ground state (S = 0) as | fla) 
(with energy E^). The block variables may be determined from these states as follows: 

U a = E^ + E<® - 2E® (2.2) 

W a -p a = E^ ) -E^\ (2.3) 

(Hi). The hopping energy between two neighboring blocks is obtained from the hopping 
energy of the neighboring sites on these two blocks. We calculate the matrix elements in 
the new states by insisting that those states are the same for both the new and the old 
Hamiltonians. There are four non-zero matrix elements for t for each bond between blocks 
a and f3: 

C = (°« h\Ht\ U Op) = (0 a h\ H t\ U Op) (2.4) 



i% = (U ip I H t | O a Up) = (U ip\H t \ O a Up) (2.5) 

t% = (U« 0/9 | St | |a T/9> = (Ua 0/9 | H t \ | Q ] p) (2.6) 
taP = (tlal/9 I H t | |a Tl/?) = (tlal/9 | #t | la tl/?)- (2.7) 

There are 3 bonds allowing hopping between any two star blocks a and (3 (also in the square 
blocks; see Fig. 1). We sum these three hops to get t^L The above 4 f Q( g's are then averaged, 
— \ J2t=i i^ap-> f° r one connection 4 between two blocks a and /3. 
(zt>). Now we average over an ensemble of random configurations to determine the 
distribution of the parameters in the block Hamiltonian. Since Q2.2|) is always positive (as 
verified numerically), we simply use its mean (denoted as U) as the renormalized on-site 
repulsion between the block particles. The renormalized chemical potential jl is defined by 
the mean of the right hand side of (|2.3|) . This implies that the renormalized random potential 
{lV a } has zero mean. However, the distribution of the block potential is in general different 
from the original one. Here we choose to keep track of only the first two moments, the mean 
(= 0) and the variance (= W); we thus map the renormalized distribution of the {W{\ back 
onto an (independent) gaussian form. 

The determination of the block hopping parameter is more subtle. Consider a block (of 
any size) without disorder. In general, the ground state is degenerate, where the degeneracy 
is related to the symmetry of the lattice and to the fermi statistics of the electrons. Any 
amount of disorder breaks this lattice symmetry and therefore lifts the degeneracy. As 
(degenerate) perturbation theory shows, depending on the configuration of the random fields 
{VFj}, the sign of the t a ps can be either positive or negative. This causes frustration when 
the product of t a ps around a closed path is negative. This is different from the boson case 



considered in Ref. |TTJ where the ground state of the pure system is non-degenerate and the 
kinetic energy is unfrustrated by site-disorder. To take this effect of frustration into account, 
one has thus to keep track of the lattice structure. More specifically, for the 'star' block 
considered here, we first build a square lattice which consists of rib = 125 coupled star-blocks 



(n s xrib = 625 sites). Step (in) is then performed to obtain the corresponding block hopping 
parameters. This explicit lattice structure enables us to compute the frustration index 
(defined by the ratio of number of frustrated plaquettes to the total number of plaquettes). 
It turns out that, starting from a uniform hopping constant tij = t, the RG described 
above will randomize the hopping parameter and frustrate the kinetic energy. Regardless of 
the starting configuration, the frustration index for the block system is always near 0.5. We 
approximate the block hopping parameters with an independent gaussian distribution which 
is determined by their mean t ave and variance t var . In the actual calculation, we typically 
average over 5-10 such lattices. Due to the symmetry of the square lattice, one can always 
choose the mean of the hopping parameter t ave to be positive. Under the RG iterations for 
the 2D problem, t ave decreases rapidly, reflecting the frustration. 

(v). Using the new set of parameters one can repeat the process described above and 
study the flow of the parameters under the RG iteration. Physical phases are identified from 
the stable fixed points of the RG. 

The Hamiltonian ( |2.1| ) is characterized by four independent parameters, which may be 
chosen as t ave /W, U/W,fl/W, and t var /W. Among the others, t var /W is found to always 
renormalize to zero; hence we shall not discuss this parameter further in the following. For 
simplicity we label the relevant dimensionless parameters as follows: U /W = U /W gives the 
strength of the repulsion; the flow of t ave /W = t/W indicates insulating (— > 0) or metallic 
(— > oo) behavior; and fl/W = ^i/W gives a dimensionless measure of the chemical potential. 

The choice of the block states described in (ii) needs some more explanation. We trun- 
cate the Hilbert space of the block by choosing four low lying states such that the block 
Hamiltonian and the site Hamiltonian have the same form. However, there are different 
possibilities of choosing these 4 states. The simplest possibility (implemented in Ref. || at 
half filling) is to always pick the same N a = N for every block a. This "fixed-n" (where n 
is defined by n = N a /n s ) procedure is the most artificial of the procedures we have used; 
however, one might argue that it is adequate for incompressible states. 

A second ( "fixed-//' ) procedure is to fix the chemical potential // rather than the density 
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n. This procedure can only work however if the chosen value for fi corresponds to a density 
which can be represented on the blocks by an integer particle number; otherwise, instabilities 
occur in the RG flow ||11|| . Hence we have only used the fixed-/i procedure at half filling. 
In this case, the chemical potential is known to be precisely U/2 and can therefore be set 
at the beginning of the RG iteration. Since the distribution of the random potential {W{\ 
is symmetric with zero mean, the statistical fluctuations will then preserve the average 
density at 1 and jl at U/2 by averaging over the random configurations. Since the density 
is allowed to vary from block to block, this method allows one to study compressible as well 
as incompressible states. 

The possibility of allowing the density to fluctuate in such a real space RG scheme is 
unique to disordered systems. In a pure system, all blocks are identical and N a will thus be 
chosen the same for all blocks. The particle number fluctuation in any given region (of any 
size) is thus one, so that the RG can only describe an incompressible state ||1 1|| . Thus one 



expects that such a real space RG scheme for pure systems will be mostly applicable at half 
filling, where it is indeed incompressible [[2tJ. 

We have used a third procedure to study the physics in a region around half filling, 
by allowing the chemical potential to also flow in the RG iterations. This enables us to 
explore the RG flow in the full 3D parameter space (t/W, U/W, (jl/W). In the absence of 
disorder, this would not be possible, since the density cannot vary in any block. In the 
presence of disorder, one selects N a which minimizes the energy in a given block a, and thus 
allows n (where again the overbar mean disorder average) to vary continuously along with 
fi. However this flow must fail at high and low densities, where the small size of our blocks 
imposes strict upper and lower bounds on the range of densities which can be handled by 
the method — for example, the density in a star block cannot fall below 0.2 nor exceed 1.8. 
Hence with this procedure it is necessary to follow the flow of n as well, and to discard 
flows when n saturates at its upper or lower bound. With this procedure it is sometimes 
convenient to parameterize the chemical potential as jx/U = ii/U since it is in terms of this 
parameter that we can locate half filling (n = 1 at \i/U = 1/2). 
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III. RESULTS 



We start by considering the Anderson model ({7 = 0). In this case, with either the fixed n 
or the fixed « method, the only meaningful parameter of the system is t/W . Either procedure 
gives the same qualitative result: we find two fixed points, at t/W = and at t/W = oo. 
The fixed point describing the pure system (t/W = oo) is unstable, with flow towards the 
attractive (insulating, t/W = 0) fixed point. Thus we find that, for noninteracting fermions, 
disorder is always relevant, in agreement with the prediction from scaling theory [rB,rB|. 

Next we consider the interacting case {7^0. In Fig. 2, we show RG flow diagrams in 
the two-dimensional parameter space (t/W,U/W). Fig. 2a is obtained at half filling. Here 
we again find that the two RG procedures (fixed-n and fixed-//) give similar results. Fig. 2b 
depicts RG flow at densities away from half filling, obtained using the fixed-n procedure. 
The possible densities are n = 1/5,3/5 for the star blocks, and n = 1/9,3/9,5/9,7/9 
for the 3x3 square blocks (apart from those which may be obtained using particle-hole 
symmetry). They all give qualitatively the same flow diagrams. At half filling, Fig. 2a 
shows that, apart from the unstable fixed points describing the noninteracting (U = 0) and 
pure (t/W = oo) phases, there are two stable fixed points at (t/W = 0,U/W = oo) and 
at (t/W = 0, U/W ~ 1.3). Between these two phases is a separatrix which terminates at a 
repulsive fixed point (t/W = 0, U/W = (U/W)* m 7.3). Away from half filling and at fixed 
density, the RG flow has only one stable fixed point, at finite U/W. We note that, with 
respect to the noninteracting system, U is relevant at all the fillings we examined. 

We next ask, what is the nature of the various phases revealed by the stable fixed points 
in Fig. 2? Since t/W renormalizes to zero in all the cases, there is no metallic state. However, 
the nature of the insulating states needs some elaboration. We consider first the half-filled 
case. For the pure system (Hubbard model) at half filling, it is believed that the system is 
always insulating with long-ranged antiferromagnetic order for any finite U > 0. However, 
while the disordered system is always insulating, the physics responsible for the insulating 
behavior may vary. This situation is best illustrated in the essentially exact calculations on 
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the infinite dimensional Hubbard model ||21|| : at low but above the Neel temperature, there 
is a critical value of U — U c , beyond which a gap opens up in the quasi-particle spectrum and 
the system changes from a metal to a Mott-insulating state. Such a paramagnetic solution 
persists down to T = 0, although it becomes unstable at low temperature and the true 
ground state is antiferromagnetic for any finite value of U. Thus one may expect that, while 
for small U the insulating state is a result of the delicate (antiferromagnetic) correlations, 
at large U it is simply due to the large energy cost for double occupancy Upon introducing 
disorder, the Mott transition (masked by the antiferromagnetic long range order in the pure 
case) is revealed, but the corresponding metallic state now becomes insulating also. Hence 
we interpret the state described by the fixed point at {U/W) — > oo as the Mott-Hubbard 
phase, while that associated with the fixed point at finite U/W we will call the 'Anderson' 
phase, since it is expected (to be corroborated below) to be gapless. 

Remarkably, this phase diagram is quite similar to that for the ID case, which was 
calculated at half filling in Ref. || using the fixed-n method. (We have obtained the same 
picture for the ID Anderson-Hubbard problem using the fixed-/! method at half filling.) 
The value of the unstable fixed point separating the two insulators is about (t/W, U/W)* = 
(0, 7.3). This value for the 2D system is very close to the fixed point obtained by Ma for 
the ID case [(t/W,U/W) = (0,8.3)], and to the slope of critical line for the opening of a 
compressibility gap in the (U c , W) plane obtained by Dommguez and Wiecko (DW) [19| 
for the 3D case, U c = 6.7W(W/t — > oo). As mentioned before, in getting the hopping 
parameter tf, between two blocks, we take an arithmetic average over 4 t a ps to take account 
of frustration. The t a ps tend to be of varying sign (for reasons discussed earlier) and hence 
to cancel each other when we average, so that the parameter t/W rapidly approaches zero 
as the RG iteration proceeds. The slope of the separatrix between the two insulators is, 
therefore, nearly zero (unlike the ID case ||). 

The stable fixed point found (with fixed density) away from half filling can be interpreted 
naturally as a fixed point describing the Anderson insulating phase. Although this phase 
diagram is obtained using the RG procedure for fixed n — which is more appropriate for 
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incompressible states — we believe that this Anderson phase is actually compressible from 
the physical point of view. The compressibility of this phase at or near half filling may be 
investigated through the RG flow in the 3D parameter space (t/W, U/W, fi/W). 

We have studied the compressibility in our numerical RG calculations using two methods. 
One, which is purely heuristic, is to stop the calculation after a single iteration, associating 
the renormalized values of density and chemical potential with the (fixed) values of U, t, 
and W which were input. This method — which can not probe the long wavelength physics 
seen from repeated RG iterations — nevertheless gives surprisingly good results, possibly due 
to the large size (625 sites) of the finite lattice which we used, coupled with the further 
averaging over disorder. Results obtained using this method are plotted in Fig. 3. We 
see that (thanks to the disorder and the averaging) the density n flows smoothly with the 
chemical potential, with two significant exceptions. One exception occurs when the density 
saturates at the maximum or minimum value allowed by the finite size of our block (i.e., for 
the star block, 0.2 and 1.8). This saturation marks a limit beyond which our method gives 
meaningless results. The other departure from smoothness occurs at n — 1, for sufficiently 
large U, and is due to the abovementioned incompressibility. The incompressibility (the 
Hubbard gap) is broadened with increasing Coulomb repulsion (Fig. 3). In the inset we 
plot the Hubbard gap Afi/t as a function of U/t, with fixed t and W. We can see that the 
Hubbard gap increases linearly with U. The slope a [= d<y ^ ) of the three curves is about 
a ~ 1.0, which is consistent with what is expected, and also with the result obtained by 
DW p|. 

From the inset of Fig. 3 we can write Afi = a(U — U c ) at fixed W. Then in (/j,/U,U/W) 
parameter space the gap can be described (for fixed t) by 

^ = fOT (U>U ). (3.1) 

Taking the large U limit, the gap Afi/U will approach a in fi/U and U/W space. Hence 
our second method for studying the behavior of the incompressibility is to follow the RG 
flow in the 3D parameter space {t/W, U/W, fi/U), distinguishing however those flows which 
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remain pinned at n = 1 from those which do not. Fig. 4 is a 3D flow diagram inside the 
limits at which the density saturates, but projected out onto the 2D plane (U/W, li/U). 
(We project out the behavior of t/W since it is the most predictable, always flowing to 
zero.) In this figure we see that the Mott-Hubbard phase (shown by the shaded area) with 
n = 1 is bounded by equation (|3.1|). We can also see the two fixed points in the plot, which 
correspond to those in Fig. 2a. For U/W — > oo we have a "fixed bar" (i.e., a line of fixed 
points in fi-t-U space) which attracts the flows inside the Mott-Hubbard phase. (This fixed 
bar is of course a fixed point in n—t—U space, with n = 1.) In the Anderson phase we can 
see that the parameter y,/U flows to the value for half filling li/U = 1/2. Hence we find 
that there is a finite region in density, around n = 1, in which the system is a compressible 
(Anderson) insulating phase and is characterized by the fixed point at half filling. 

Since the discovery of high T c superconductivity, there have been suggestions that the 
dimensionality alone will invalidate the Fermi liquid theory and make a 2D interacting system 



a highly correlated one [23], as it does for its ID counterpart. While for the 2D Hubbard 
model near half filling there is little doubt that the system is indeed highly correlated, 
controlled perturbative expansions suggest that at low density a system of interacting 2D 
fermions can be well described by the conventional Fermi-liquid theory [[23]]. If this is the 
case, the insulating state at low density would be simply due to the localization of the quasi- 
particles. From this point of view, one would expect that the (compressible) insulating 
state near half filling will be quite different from that at low (or high) densities, so that 
there should be an additional stable fixed point describing such a 'conventional' Anderson 
insulating state, distinct from the 'highly correlated' Anderson insulating state controlled by 
the fixed point at half filling. Hereafter we shall call the 'conventional' Anderson phase an 
'Anderson- Fermi' insulator; the (presumably) highly correlated Anderson phase described 
by the Anderson fixed point at half filling we call the 'Anderson-Luttinger' insulator. We 
choose the latter name since our RG study shows a strong resemblance between the ID 
and 2D systems near half filling, and since generic (pure) ID systems are described by the 
highly correlated 'Luttinger liquid' f22 |. In our studies we have found no evidence for an 



13 



Anderson- Fermi insulating phase characterized by a high- or low-density fixed point. We 
note however that we cannot rule out the existence of such fixed points, since our method 
is only reliable in a finite range of densities around half filling, and so may be incapable of 
detecting these uncorrelated phases. 

In Fig. 5 and Fig. 6 we show the full flow diagrams in the 3D parameter space. Taking 
account of particle- hole symmetry, we omit the region n < 1. We also omit U < (which 
is expected to give different physics) and t < which is trivially related to t > 0; hence 
we are left with one octant of the full space. In each plane the flows are the projection 
of 3D parameter flows. The incompressible Mott-Hubbard phase is lightly shaded in the 
(n/U, U/W) plane. Our RG approach fails in a high- and a low-density region; the former 
is shown in Figs. 5 and 6 with dark shading. In this dark region we can imagine either two 
fixed points (one in the high-density region and the other one in the low-density area), or 
none. In the case of no fixed points in the dark area (Fig. 5) — or if the noninteracting fixed 
points are unstable to any finite U — the RG parameters outside the Mott-Hubbard phase 
flow to the one fixed point, so that there is only one ("correlated") Anderson phase in the 
Anderson-Hubbard model. On the other hand, if there are two stable fixed points (at high 
and low density) (Fig. 6), the Anderson-Hubbard problem will have two different Anderson 
phases as discussed above. We include both Figs. 5 and 6 because we do not believe that 
we can distinguish these two scenarios within the limitations of our method. 

IV. SUMMARY AND DISCUSSION 

Using a quantum real space renormalization group method, we have obtained the phase 
diagram of the 2D Anderson-Hubbard model near half filling at T = 0. A test of our method 
for noninteracting fermions shows the instability of the metallic phase for any nonzero dis- 
order W, in agreement with commonly accepted results and hence providing some evidence 
that our method is qualitatively reliable. By studying the renormalization of chemical 
potential and the corresponding flow of particle density, we were able to estimate the com- 
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pressibility gap A/i in the Mott-Hubbard phase as a function of U and W. We found that 
A/i increases proportional to the interaction U (at constant W) above a critical value (U/t) c , 
with the constant of proportionality a about 1, and that A/i decreases with increasing W 
(at fixed [/). Our results here are also in good agreement with those obtained by other 
methods, which gives us further confidence in our approach. 

Our studies show that there is no metallic phase for the 2D Anderson-Hubbard model, 
for any finite value of the random potential W and of the repulsive interaction U between 
the electrons. The interplay between interaction and disorder yields two insulating phase 
at half filling: an incompressible Mott-Hubbard insulator and a gapless Anderson insulator. 
The phase diagram strongly resembles that for the corresponding ID system. Away from 
half filling, the insulating phase is always gapless, and its properties are controlled by the 
fixed point describing the Anderson insulator at half filling. We characterize such a highly 
correlated insulating phase as the 'Anderson-Luttinger' insulator. U is relevant with respect 
to the noninteracting fixed point for all the cases we have considered. We would like to 
emphasize, however, that the relevance of U itself does not constitute evidence for the 
existence of the 'Anderson-Luttinger' insulator. Since the noninteracting disordered system 
is described by the localized (t/W = 0) fixed point, one may expect that, at least in the low 
doping (near half filling) case, any interaction will be relevant regardless of the properties of 
the corresponding pure system, although there is a possibility that short-ranged interactions 
such as the on-site U studied here are not relevant in the dilute (low or high density) 
limit. The dilute fixed point, which presumably describes the conventional 'Anderson- Fermi' 
insulator whose physics is related to that of localized noninteracting fermions has not been 
found within our approach, which is however limited to a range of densities around half 
filling. 

The picture that we find at half filling is not unexpected: the instability of the nonin- 
teracting fixed point, the consequent finite-?/ (gapless) fixed point, and the opening up of 
a gap at U c with flow towards U — oo in the Mott-Hubbard phase. The one feature of 
our results that is perhaps somewhat unexpected is the existence of a finite region in n (or 
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/i) around half filling, which is dominated by the n = 1 fixed point. One could imagine a 
different result, namely, that, like the Mott-Hubbard phase, the 'Anderson-Luttinger' phase 
is well-defined only at or close to half filling, becoming unstable as n deviates from this 
region and flowing towards a 'dilute' fixed point. In other words, one could imagine that 
5fi = fi — U/2 is relevant when its magnitude becomes nonzero or sufficiently large, which 
is not what we have found for the region allowed by our RG scheme. 

It is conceivable that this result might be due to an artifact of our method. That is, one 
might conjecture that the flow towards half-filling reflects only the stability of the algorithm 
at half filling, rather than the stability of the thermodynamic phase. Although we see no 
reason for this to be the case, we cannot rule out this possibility. We do however gain 
some confidence in the results of our RG method, in the case where we allow the chemical 
potential to flow, from the good agreement of these results for the Hubbard phase with 
existing results obtained by other methods. We note that these results (Figs. 3 and 4) were 
all obtained using this algorithm. 

We therefore assume that the correlated 'Anderson-Luttinger' insulator indicated by 
our results is in fact the true ground state of the 2D Anderson-Hubbard problem for some 
region around half filling. This suggests a number of directions for future work. It is clearly 
important to try to clarify the nature of this phase, in both the 2D and the ID problems, 
for instance by calculating density-density or magnetic correlation functions and their RG 
flow. Furthermore, if indeed the physics for this disordered problem around half filling is 
described by the Anderson fixed point at half filling, one can expect to gain significant 
information about the lightly doped case by directly studying the half-filled case (where, 
for instance, there is no "sign problem" in quantum Monte Carlo simulations). It would 
also be of considerable interest to extend this work to the 3D problem, where one expects 
a metallic phase, and metal-insulator transitions of various types ||. Unfortunately, the 
smallest isotropic 3D block with an odd site number (3x3x3) [|17]] is far too large for exact 
diagonalization. Since the 3-(spatial)-D problem is of interest both in its own right, and as 
a further test of the present 2D results, we believe that the problem of extending our real 
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space RG technique for disordered systems to the 3D case merits some further effort. 
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FIGURES 

FIG. 1. Two types of blocks used in our RG approach for the square lattice, (a) star lattice 
(b) 3 x 3 lattice. 

FIG. 2. Flow diagrams of the 2D Anderson-Hubbard model as obtained by our RG approach, 
at fixed filling of the lattice, (a) For half filling we see two insulating phases: a Mott-Hubbard 
phase at large U /W and an Anderson (gapless) phase at smaller U. (b) Away from half filling we 
see only the Anderson phase. 

FIG. 3. Density versus fi/U with different U/W for fixed t(= 0.5) and W(= 1.0). Each curve is 
shifted to the right, for visual purposes, by 0.5 as U/W is increased by 1. Density is pinned at n = 1 
over a finite range in \x as U is increased, (inset) The incompressibility A/j, in the Mott-Hubbard 
phase is plotted as a function of U in units of t, for fixed W and t. For comparison, we include 
some results for the ID lattice. 

FIG. 4. Projection on a plane of flows of 3D RG parameters {t/W, U/W, fi/U}. In our projec- 
tion the flow of t/W is not shown, since this parameter always flows to zero. The Mott-Hubbard 
phase (in which the density is pinned at 1) is shaded. The unstable fixed point (U/W = 7.3, 
fi/U = 0.5) marking the boundary of the Mott-Hubbard phase is marked with a small circle; the 
stable one is located at (U/W,n/U) = (1.3,0.5). 

FIG. 5. 3D parameter flow diagram. In the dark area, where the density is saturated and so 
our method gives no information, we are assuming the chemical potential flows toward half filling. 
The Mott-Hubbard phase is marked with light shading. 

FIG. 6. 3D parameter flow diagram, but with an alternative hypothetical scenario from that 
shown in Fig. 5. Here we assume the existence of attractive fixed points in the dark areas (high 
and low density), indicating the presence of "uncorrected" insulating phases distinct from the 
"correlated" phase we find around half filling. Since our method fails in the dark shaded area, it 
cannot distinguish between the picture shown here and that in Fig. 5. 



20 



